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ABSTRACT 

NGC 326 is one of the most prominent X- or Z-shaped radio galaxies (XRGs/ZRGs) and has been 
the subject of several studies attempting to explain its morphology through either fluid motions or 
reorientation of the jet axis. We examine a 100 ks Chandra X-ray Observatory exposure and find 
several features associated with the radio galaxy: a high-temperature front that may indicate a shock, 
high-temperature knots around the rim of the radio emission, and a cavity associated with the eastern 
wing of the radio galaxy. A reasonable interpretation of these features in light of the radio data allows 
us to reconstruct the history of the AGN outbursts. The active outburst was likely once a powerful 
radio source which has since decayed, and circumstantial evidence favors reorientation as the means 
to produce the wings. Because of the obvious interaction between the radio galaxy and the ICM and 
the wide separation between the active lobes and wings, we conclude that XRGs are excellent sources 
in which to study AGN feedback in galaxy groups by measuring the heating rates associated with 
both active and passive heating mechanisms. 

Subject headings: galaxies: active — galaxies: individual (NGC 326) 



1. INTRODUCTION 

A small fraction of double-lobed extragalactic radio 
sources have long, centro-symmetric tails of radio emis- 
sion extending from the nucleus or primary lobes in 
a very different direction, thereby producing "X" or 
"Z" -shaped radio galaxies (XRGs/ZRGs). The tails, or 
"wings" , are faint, do not currently harbor a relativistic 
jet, and can be substantially longer (in projection) than 
the primary lobes. They have elicited considerable in - 
terest due to the suggestion by 'Mcrritt fc EkersI ()2002[ ) 
that they result from a recent supermassive black hole 
(SMBH) merger that reoriented the jet, but their ori- 
gin is uncertain. It has also been suggested that the 
wings result from the deflection of lobe plasma by pres- 
sure fronts in the intergalactic or intracl uster medium 
jlGM/ICM) (Leahv & Williams 1984: iWorrall et al. 
1995; Capctti ct al. 2002; Gopal -Krishna et al.ll2003l: lZiei 
2005; .Hodges-Kluck et al.n2010i) . Neither of these hy- 



potheses have yet been proven, even in a single source. 

B2 0055-^26 is a dramatic XRG/ZRG, with wings 
longer (in projection) than the primary lobes (Figure [1]). 
The wings are highly coUimated over '--^lOO kpc. Its 
host gala xy is the northern component of the dumbbell 
NGC 326 (jColla et al.lll975[ ) and the ent ire system resides 
in th e galaxy cluster Zw 0056.9-1-2636 (jZwickv fc Kowall 
Il968f ). Hereafter, we adopt the convention in the litera- 
ture and refer to the whole system as NGC 326. The clus- 
ter at mosphere has been characterized bv IWorrall et al.l 
(|1995L hereafter, W95) on large scales using ROSAT, and 
has an average temperature of kT ~ 2 keV. The cluster 
atmosphere is substantially asymmetric (Figure [2]), sug- 
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gesting that it is not in hydrostatic equilibrium and may 
be made up of two or more merging subclusters or that 
it is a composite of two or three poorer clusters seen in 
projection (W95). A detailed multifrequency study of 
the radio emission measured the spectral index a as a 
function of position in the lobes and wings, finding spec- 
tral stee pening in the wings and away from the ends of 
the jets ijMurgia et alll2001l hereafter MOl). 

Owing to its dramatic structure (Figure [1]), NGC 326 
is frequently used as the prototypical XRG/ZRG in for- 
mation models purporting to explain the X-shaped mor- 
phology. lEkers et al.l ()1978t ). who discovered the wings 
but not their true length, attributed them to ongoing 
regular precession o f the relativistic jet infi ating the pri- 
mary lobes. Later. [Merritt fc Ekerd ()2002i) argued that 
the striking X-shaped structure necessitated rapid secu- 
lar precession. Noting that black holes are difficult to re- 
orient via external forces other than accretion disk insta- 
bilities, IMerritt fc EkersI (|2002t ) suggested instead that 
the X-shaped structure is the result of the final coales- 
cence of a SMBH binary. Supposing that the components 
of the binary have misaligned angular momenta, the fi- 
nal spin axi s (and jet axis) will be nearly instantaneously 
reoriented. iWirth et al.l (Il982| ) also argued in favor of a 
single rapid jet reorientation caused by a close encounter 
between the constituents of the dumbbell galaxy. 

Based on the large-scale, hot, asymmetric cluster at- 
mosphere, W95 found that the buoyancy of the lobe ma- 
terial could explain the wings provided that they were 
driven by rapid backflows expected near the terminal 
shocks of the radio jet heads. This scenario is funda- 
mentally different from those above in that it presumes 
a st able jet axis. In support o f the backflow hypothesis, 
we (|Hodges-Kluck et "al1l2010D reported a local asymmet- 
ric enhancement of the hot gas in a preliminary analysis 
of an archival 100 ks Chandra X-ray Observatory (CXO) 
image on the scale of ^50-75 kpc where the radio beam 
is aligned with the major axis and the wings with the 
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Fig. 1. — Left: 1.46 GHz VLA C Array image. The contour levels are 1, 2, 4, 8, 16, 32, 64, and 128 mjy/beam. The restoring beamsize 
is 16. "7 X 14. "4. Right: 1.55 GHz VLA A Array. The contour levels are 0.1, 0.2, 0.4, 0.8, 1.6, and 3.2 mJy/beam. The restoring beamsize 
is l."35 X l."15. Note that the images are on different size scales, and the full extent of the wings is not seen at right. 



m inor axis. This is als o consistent with the expectation 
of lCapetti et all (|2002[ ) that X-shaped sources are aware 
of the geometry of their environments. However, it is not 
clear how the wings would grow so long in these models. 

In this paper, we carefully consider the high-resolution 
CXO data in order to shed light on the history of 
NGC 326, evalauting possible hydrodynamic hypothe- 
ses and characterizing the interaction between the radio 
galaxy and its environment, which could not be resolved 
in the W95 study. NGC 326 is located at a redshift z = 
0.0474 where 1" = 0.92 kpc in the Wilkinson Microwave 
Anisotropy Probe cosmolog y (Hp = 71 km s~^ Mpc~^, 
Qa = 0.73, and = 0.27. ISpergel et al.|[2007l ). We use 
a Galactic photoelectric absorpt ion in all spectral fitting 
of TVh = 5.69 X 10^0 cm-2 (Kal berla et al.l200 5). The to- 
tal extent of the cluster atmosphere in the ROSAT image 
is about 20' on a side (1.1 Mpc). 

The remainder of the paper is organized as follows: in 
Section 2 we describe the CXO observation, data, and ba- 
sic features, then in Section 3 we use the X-ray and radio 
data to derive quantities of interest related to the radio 
galaxy. In Section 4 we use these data to reconstruct 
the history of the radio galaxy and argue that NGC 326, 
as well as other XRGs, are interesting sources with re- 
gard to feedback from active galactic nuclei (AGNs) . We 
summarize our results in Section 5. 

2. OBSERVATIONS 

NGC 326 was observed on 2006-09-02 for 94.4 ks by 
CXO using the Advanced CCD Imaging Spectrometeo 
(ACIS) with the dumbbell galaxies located at the nom- 
inal aimpoint on the S3 chip (obs. ID 6830). We previ- 
ously pu blished these dating look ing only at gross mor- 
phology ( Hodges-Kluck et al.l [2010). The data were re- 
duced with the Chandra Interactive Analysis of Obser- 
vations (CIAO v4.0) software, starting from the level=l 

* See |http://cxc.harvard.edu/proposer/POG/pdf/ACIS.pdf| 



file rather than using the pipeline level=2 file. No back- 
ground flares in the 0.3 — 10 keV bandpass were detected 
in the lightcurve of low surface brightness regions on the 
S3 chip or empty regions on other chips (the entire S3 
chip is covered by some cluster emission). Fo r spectral 
fitting we used XSPEC vl2.7.0 (IArnaudlll996D. w ith the 
apec model as our thermal model ( Smith et al.l l200"ll) . 
Images produced from the reduced events file were ex- 
posure corrected using an exposure map generated in 
CIAO, and questionable pixels at the very edge of the 
chip which result from this process are excised from our 
analysis. For analysis of the hot cluster atmosphere, 
point sources are excised from the image using the CIAO 
wavdetect tool. We also exclude the resolved interstellar 
media (ISM) of the dumbbell galaxies. Since the cluster 
emission extends north onto the S2 chip, we follow an 
analogous procedure for these data. 

In addition to the CXO data, we used archival NRAC0 
Very Large Array (VLA) radio maps presented in MOl. 
These include maps at 1.4, 1.6, 4.8, 8.5, and 14.9 GHz 
produced with the NRAO AIPS software. Re-reduction 
of one (post-calibration) dataset with the newer CASA 
software produced no significant differences. As de- 
scribed in MOl, the angular size of the source is com- 
parable to the size of a VLA antenna, so two point- 
ings are combined to form the final image. The dumb- 
bell system NGC 326 itself has been observed with the 
Hubble Space Telesc ope (HST) as a snapshot observation 
(jCapetti et al.ll2000H . and there is archival ROSAT data 
(W95). 

In the ROSAT PSPC image (Figure [2]) the X-ray peak 
coincides with the radio galaxy but is not near the center 
of the cluster (W95). The ROSAT image does not have 
point sources subtracted; the bright spot near the south- 
east edge of the CXO S3 chip is a point source subtracted 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Fig. 2. — Left: ROSAT 0.5-2.4 keV image smoothed with a cr = 3 Gaussian kernel (no point source subtraction). The black boxes show 
the outlines of the CXO chips with cluster emission. The radio contours at 4.9 GHz are overlaid in blue. Right: CXO 0.5-5 keV image of 
the S3 chip with point sources subtracted and binned to 16 X 16 pixels, then smoothed with a cr = 3 Gaussian kernel. The radio source is 
shown in blue as on the left. 



in the CXO image shown in Figure [H The CXO image 
reveals a significant local anisotropic enhancement in the 
soft X-rays near th is peak elongated in the d irection of 
the primary lobes (|Hodges-Kluck et al.|[20Tol ) that can- 
not be seen in the ROSAT image. In fact, the ROSAT 
X-ray peak is cospatial with the northern lobe of the 
radio galaxy and not with the dumbbell galaxies; these 
galaxies are visible in the CXO image but on scales much 
smaller than the ROSAT point spread function (PSF). 

2.1. The ISM of NGC 326 

The CXO image resolves the constituents of the dumb- 
bell NGC 326 (Figure m smoothed X-ray contours over- 
laid in green). The northern component is the host 
galaxy of the radio source and has an obvious strong 
point source at its center. As in the HST images, 
the northern comp onent is the brighter of the pair 
(jCapetti et al.l [20021 ). Spectral fits for the two galaxies 
are described below. 

Northern Component — In iHodges-Kluck et al.l ()2010f ) 
we measured a temperature of kT ~ 0.7 keV in this 
galaxy both including {kT = 0.68 ± 0.06 keV) and ex- 
cluding {kT — 0.7 ± 0.1 keV) the central point source, 
which we fit using a power law of F = 1.3 ± 0.4. Metal- 
licity was frozen at solar abundances. This fit produces 
a = 11-1 for 16 degrees of freedom (20 counts bin~^). 
The (unabsorbed) model luminosity of the thermal com- 
ponent is Lx ~ 7.4 X 10^° erg s~^, whereas the X-ray 
luminosity of the power law is Lx ~ 5.4 x 10"'" erg s~^. 

An acceptable fit can be obtained using a 2-T fit with 
kTi = 0.7 keV and kT2 > 3 keV. Presumably a higher 
temperature contribution would come from the projected 
ICM, but this is not a satisfactory explanation because 
the excess high-energy photons that require the higher 



temperature or power law component are concentrated 
inside the PSF at the location of the radio point source. 
Meanwhile, photons outside the half-power distance of 
the PSF are preferentially softer, suggesting that they 
are really part of the extended ISM. An AGN origin for 




0.25" 0" -0.25" 

Fig. 3.— HST WFPC2 image of NGC 326 (median-filtered to 
remove cosmic rays) overlaid with X-ray contours (green) and 
1.55 GHz contours (black, primarily showing the radio jet). The 
ISM of both components of NGC 326 are visible and resolved by 
CXO. The northern component is brighter in both the optical and 
X-ray, and is the host of B2 0055+26. 
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the harder X-ray emission is therefore favored. 

Southern Component — The southern component is 
fainter and has no obvious nuclear point source. Only 
154 counts are detected, with a total luminosity of Lx ~ 
6 X 10^° erg s~^. Despite the small number of counts, 
the X-ray spectrum is poorly fit by a single 1-T model 
(X^ ~ 12 for 8 degrees of freedom) but a good fit is 
achieved by including a power law component with a 
frozen photon index F = 2 (x^ = 7.3 for 8 degrees of free- 
dom). There are insufficient counts to constrain F, but a 
power law interpretation is preferred as in the northern 
component since the harder photons are concentrated to- 
ward the center of the detect cell. Some power law emis- 
sion is expected from unresolved X-ray binaries (see, e.g. 
iDiehl fc Statleii[2007D . but we cannot rule out a weak nu- 
clear source due to the 1 GHz emission from the nucleus. 

A common envelope of diffuse optical emission sur- 
rounds the components of the dumbbell galaxy, but no 
similar envelope is detected in the X-rays. The radio im- 
age leaves open the possibility of interaction between the 
jet and the ISM of the southern component, but no X-ray 
jet is detected and the signal-to-noise ratio (S/N) in the 
southern component is too low to search for disturbances. 

2.2. Cluster Atmosphere 

The cluster X-ray emission is diffuse and extensive, 
covering the ACIS-S2 and S3 chips as well as some of the 
12 and 13 chips. No other chips were active for this ex- 
posure. The aimpoint for the observation was on the S3 
chip, so it is difficult to identify and excise point sources 
on the 12 and 13 chips. While bright point sources are 
still detected with algorithms such as wavdetect, faint 
point sources are distributed over many pixels and are 
difficult to distinguish from the cluster atmosphere or 
the background. Thus, we focus our efforts on the S2 
and S3 chips (the radio galaxy is wholly contained in the 
latter). Quantities are listed in Table [T] 

Large-scale structure — We generated a radial profile 
from the point-source subtracted images, centered on 
the X-ray peak near NGC 326 (Figure |4]) , subtracting 
local background from an annulus outside the fit re- 
gion. We fit a ID /3 model (a Lorentzian of form /(r) = 
Ao[l + (r/ro)^]'~^^'''^/^\ where Aq is the peak amplitude 
in counts per square arcsecond and rp is the core radius 
in arcseconds) and find a good fit for (3 — 0.53 ±0.04 and 
ro = 110 ± 10 arcsec (100 kpc). This gives a = H-IO 
for 11 degrees of freedom (using CIAO's Sherpa with 
the Gehrels variance function). Our value of /3 is con- 
sistent with the W95 value, but they find a core radius 
ro ^ 180 kpc in WMAP cosmology. This may be due to 
the larger field of view of the ROSAT image and the fact 
that our background annulus encompasses some cluster 
emission. There is also obvious anisotropy in the cluster 
on the scale of the S3 chip (Figure [5]), so our fit is only 
a good description of the atmosphere in the vicinity of 
NGC 326. Since the X-ray emission may be a superpo- 
sition of multiple clusters or groups (W95), we may be 
isolating group-scale gas associated with NGC 326. 

W95 measured a cluster temperature of approximately 
2 keV and a metallicity of 0.3 times solar from the 
ROSAT data. The CXO data is largely in agreement, 
with temperatures ranging from kT > 2 keV near the 



core to kT ~ 1.5 keV near the edge of the cluster. How- 
ever, a temperature map made using weighted Voronoi 
tessellations to adaptively bin the images into bins with 
S/N = 35 (IDiehl fc StatleH [20061 : ICaPDellari fc Copinl 
120031 ) shows a "front" of high-temperature plasma to the 
southeast of the radio galaxy (Figure |6]). The temper- 
atures in bins in this front range from 3 — 6 keV for a 
fixed metallicity oi Z = 0.3^0 (W95). The temperature 
is not well constrained, but no acceptable fits were found 
below kT ^ 2.5 keV for different bin sizes, metallicities, 
or 2-T fits. 

Another way to estimate temperatures in the front is to 
"deproject" the emission by assuming a spherical atmo- 
sphere (within 100 kpc of NGC 326) and using annuli to 
break the emission into multiple bins. In this assumed 
geometry, a line of sight through the center sees a su- 
perposition of all the cluster emission, whereas a line of 
sight near the edge sees only the emission from the out- 
skirts. Thus, each annulus is corrected for the annuli 
exterior to it, hopefully mitigating the effects of projec- 
tion. Since the high-temperature front is about 70 arcsec 
from the core and the core radius is — 110 arcsec, we 
would expect significant contamination and that the de- 
projected temperature would be higher. This is indeed 
what we find (Figure |3]), but we lack the counts to deci- 
sively bound the temperature. Since a power law model 
is not a good fit, the detection of hot plasma southeast 
of the radio source seems secure. We discuss the nature 
of this front in Section 3. 

In addition to the high-temperature front and surface 
brightness enhancement around the radio galaxy, the 
CXO map reveals a few small-scale structures unresolved 
by ROSAT. 

High-temperature knots — First, separating the data into 
soft (0.3 — 3 keV) and hard (3 — 7) keV bands reveals the 
presence of bright knots around the radio galaxy in the 
hard X-ray band (Figure [S]) . Above 7 keV, the image is 
very noisy. The presence of bright knots around the radio 
emission suggests interaction between the radio source 
and the ICM. This inference of association is strength- 
ened by the observation that the western wing appears 
to pass through a narrow channel devoid of hard X-ray 
emission in Figure [SJ Although higher signal is required 
to establish their morphology, the data suggest that these 
knots are part of high-temperature rims bounding the ra- 
dio source: the hard X-ray emission appears to avoid the 
radio emission. The presence of the brighter knots is con- 
sistent with t he detection of X-ray emitti ng rims around 
radio sources ([McNamara fc Nulsenll2007l and references 
therein). However, usually such rims are cool, whereas 
these features are distinct only in the hard 3 — 7 keV 
band. 

There are ~500 total events in the 3 — 7 keV band in 
these features. When spectra are extracted from these re- 
gions, the spectra are fit well by 1-T thermal models with 
kT ^ 3.2 keV. Given that we view the cluster in projec- 
tion and the surrounding material has kT ^ 2 keV, this is 
a lower bound on the temperature in the knots. Despite 
the small number of counts, a nonthermal origin for the 
hard emission is disfavored. If we fit the spectrum with 
a thermal model frozen at the kT ~ 2 keV temperature 
of the background plasma and add a power law compo- 
nent, we require F = 1.3^q'^ for a good fit, correspond- 
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Fig. 4. — Left: ID radial profile extracted from the ACIS-S3 chip centered on NGC 326. Point sources have been extracted as well as 
extended ISM emission. The best-fit ID /3 model is overplotted. Dashed lines denote the region covered by the high temperature bin in 
the right panel. Right: Deprojccted temperature profile for annuli 40" in width, centered on NGC 326. Temperatures and errors obtained 
using steppar with projct (phabs*apec) . The temperature in the bin associated with the high-temperature front is not well constrained, 
but must be above ~ 5.5 keV. 



ing to a spectral index a = 0.3. There are two plausible 
possibilities for nonthermal emission resulting from in- 
teraction of the ICM with the edge of a radio galaxy: 
inverse Compton emission and Fermi acceleration from 
shocks. In the latter case, we would expect a = 0.7, so 
we can rule this out. Inverse Compton emission is also 
unlikely in this case. The spectral index of the radio 
galaxy is a ~ 0.7 from 25 MHz through 1 GHz, so if it is 
synchrotron self-Compton, the boosted electrons would 
have to come from a population below 25 MHz. The 
average frequency boost is Vohs ~ I'^oT^j so this would 
imply 7 > 2 X 10^! Moreover, the knots are outside the 
radio emission at GHz frequencies, and we would expect 
to see GHz emission from these regions as well. Sup- 
posing that the inverse Compton emission comes from a 
different source, we would still expect to see synchrotron 
emission in the GHz band from the boosted electrons if 
Sknot ~ Bcq (discussion follows in Section 3). Thus, we 
conclude that the knot emission is thermal. 

Wing Cavity — In addition to positive features, an un- 
sharp mask image of the S3 chip (Figure [T]) reveals an X- 
ray decrement associated with the east wing of NGC 326 
centered about 95-10 kpc from the host. Although there 
are other X-ray decrements in the field, this is the deep- 
est. The unsharp mask uses a radius of 50 pixels, a con- 
trast enhancement of 5%, and a threshold of counts 
to apply the unsharp mask. For viewing ease we show 
the entire chip smoothed by a Gaussian kernel with width 
a — S pixels, but the chip edges enhanced by the unsharp 
mask are unphysical and should be ignored. This decre- 
ment is also plainly visible in a residual image produced 
by subtracting off the /3-model fit above (not shown). 
The decrement is in a shallow trough in the exposure 



map, but the trough spans the length of the chip and our 
exposure corrected image demonstrates that the decre- 
ment described here is substantially deeper than the 
trough. We tentatively identify this decrement as a "cav- 
ity" associated with the eastern wing. 

Computing the significance of the decrement is not 
straightforward because of the obvious anisotropy of the 
cluster emission on scales larger than the chip. Subtract- 
ing off the /3-model reveals a gradient in the residuals 
apart from the decrement. There is no obvious way to 
account for this emission, since increasing the number 
of ad hoc model components risks overconstraining the 
problem. However, we can at least ask whether a cavity 
at the location in t he cluster could be detected by using 
the expressions in iDiehl et al.l (|2008| ). Equation 13 in 
their paper gives the S/N of a. cavity of a given size and 
at a given distance from the cluster core as a function 
of /3, ro, kT, no, and observational parameters such as 
the detector effective area and exposure time. For our 
cluster {kT - 2 keV, no - 3 x lO'^ cm'^, /3 ~ 0.5, 
ro ~ 100 kpc), a cavity at the position of the decrement 
would require a radius of Tcav ^21 kpc in order to be 
detected with a S/N > 3.0 in our CXO exposure. 

This radius is consistent with the X-ray decrement seen 
in Figure [7] and the width of the eastern wing on the sky; 
a real cavity is unlikely to be larger than this. Hence, the 
decrement seen in the unsharp mask and residual image 
is consistent with the detection of a cavity. A deeper 
exposure is required to firmly establish its presence: a 
higher effective area (i.e., XMM- Newton) or a iobs ^ 
400 ks exposure with CXO would easily detect a 20 kpc 
cavity at more than 7cr. 

3. DERIVED PROPERTIES 
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Fig. 5. — Left: 0.3-3 keV CXO image of S3 chip, binned to 16 X 16 pixels then smoothed with a cr = 3 px Gaussian kernel. 4.9 GHz 
contours are overplotted. The units of the scale are in counts and the colormap is linear. Right: 3-7 keV image processed in the same way 
as left except that the units of the scale are exponential to highlight the brightest regions in the map. 
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Fig. 6. — Temperature map made from adaptive bins with S/N = 
35 each. The bins on the S2 chip are larger owing to the lower 
surface brightness. Temperatures are measured from 1-T apec fits; 
blacked out bins represent bins for which a temperature could not 
be derived with an acceptable fit. 4.9 GHz contours of NGC 326 
are overlaid for reference. 



TABLE 1 
Derived Parameters 



Parameter 



Value 



Units 



Cluster Parameters 



ro 


110±5 


arcsec 


/3 


0.53±0.04 




Ao 


0.22±0.02 


cts/(arcsec)^ 


kT 


~2 


keV (avg) 


no 


3 X 10-3 


cm"-^ 


Po 


~ 10-" 


dyne cm"^ 


Cs 


500 


km 


Radio Galaxy 



1.1 


N. Lobe Length 


-50 


kpc 


S. Lobe Length 


-75 


kpc 




E. Wing Length 


~100 


kpc 


0.S6 


W. Wing Length 


~100 


kpc 




1 GHz tgync 


70 


Myr 







~7 





Note. — Measured and derived properties of 
NGC 326 and environs. Where errors are quoted, 
they are 90% confidence limits. Otherwise, quan- 
tities should be treated as approximate with jus- 
tification in text. The radio structures arc seen 
in projection and may be longer than the given 
values. 



Properties derived in this section are summarized in 
Tabled 

Sound Speed and Terminal Buoyant Velocity — Whereas 
powerful radio galaxies grow hypersonically as their jets 
plough through the tenuous IGM, the low-power FR I 
sources are thought to expand transonically or slower 



over the bu lk of their lifetime (for a typical model, see 
iLuo fc Sadle r 2010). The remnants of dead radio galax- 
ies are also thought to rise buoyantly in the ambient 
medium, and the growth of the wings of XRGs/ZRGs 
may be subsonic depending on the mechanism. Hence, 
the sound speed is a useful yardstick to apply to these 
sources. 
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Within 100 kpc of the cluster core, the average temper- 
ature in regions other than the high-temperature front is 
about kT « 2 keV, yielding a sound speed 



Cs = 



< 500 km = 0.5 kpc Myr' 



(1) 



in the approximation of an isothermal relaxed atmo- 
sphere and /i = lAniH- 7 is taken to be 5/3. 

The terminal buoyant velocity of a light spherical bub- 
ble suspended in a medium is given by 



^buoy 



16 Tf, 



(2) 



where Vf,, A^, and rb are the volume, cross-sectional 
area, and radius of the bubble respectively, g is the local 
gravitational acceleration, Cd ~ 0.75 is the drag coeffi- 
cient, and ro « 100 kpc is the core radius of the atmo- 
sphere. For the cavity at the edge of the eastern wing, 
r ~ 100 kpc and r?, ^ 20 kpc, so 



I'buoy ~ 0.8c, = 0.4 kpc Myr~ 



(3) 



Eqm partition estimate of B — Following iTavecchio et al.l 
(|2006f) . we estimate the equipartition field Sgq from the 
1.4 and 4.9 GHz maps of the active lobes: 



2 X 10 



-5 172a 
' mill 



(2a - l)c{a)V 



(4) 



where S is the Doppler factor, 7inin is the minimum elec- 
tron cutoff energy, is the synchrotron frequency (we 
assume the observed frequency Uq = Vg), Lg is the syn- 
chrotron lumi nosity in vo l ume V , and c(a) is a numer- 
ical factor in iPacholczvkl (|1970[ ). The lobe volume is 
approximately V ~ irr'^h, where r and h are the mea- 
sured half-width and length of the lobes respectively: 
r ~ 20 kpc and /i ~ 75 kpc (south) or 50 kpc (north). 
Hence, V ^ b y. 10^^ cm^'^. The luminosity is 



(5) 



where Dl = 210.3 Mpc is the luminosity distance. High 
quality images exist at multiple frequencies, showing no 
evidence for spectral steepening between 1.4 and 4.9 GHz 
(q^i.4-4.9 ~ 0.7; MOl). Since B^^ depends on the quan- 
tity i/gLsivs), estimates of Bcq at 1.4 and 4.9 GHz must 
come out to be the same. 

The flux at 1.4 GHz is i^s,i.4 ~ 0.86 Jy, while at 
4.9 GHz 4 9 « 0.42 Jy. These values correspond to 
Bcq(5 ^ 7/iG at 1.4 GHz and -6 ^^G at 4.9 GHz. Since 
the luminosity in these images is dominated by lobe emis- 
sion, presumably S ^ 1. Hence, we adopt 



7fiG 



(6) 



as the equipartition field in the lobes. 

From the synchrotron frequency, — 4: x lO^Bj^ Hz, 
we then find 7 ^ 7100, yielding a synchrotron cooling 
time of 



sync 



24007^:^51^ Myr 



70 Myr 



(7) 



for the 1.4 GHz electrons, where 74 is 7 in units of 10** 
and B-Q is B in units of fj,G. For the 4.9 GHz electrons, 




Fig. 7. — Unsharp mask image made using a radius of 50 pixels, 
a 5% contrast increase, and a threshold of counts. Zero is rep- 
resented by the uniform orange color outside the chip. The image 
is binned to 16 X 16 pixels and smoothed by a Gaussian kernel of 
cr = 3 px. The enhancement at the chip edges is not physical, and 
the significance of the decrement associated with the eastern wing 
is based on an unsharp mask on only those pixels making up chip 
S3. The negative associated with the wing is by far the deepest in 
the image. 



isync '~ 37 Myr. For our purposes, the synchrotron cool- 
ing time is more relevant in the wings, which will fade 
from view on this timescale without a fresh injection of 
electrons. Since the wings are bright at 4.9 GHz, the 
most recent injection must be no more than about 40 Myr 
old. 

High-temperature Front/Bow Shock — The high- 
temperature front visible in the temperature map 
(Figure ^ and evident in the deprojected temperature 
profile (Figure |4]) follows the contours of the soft X-ray 
emission to the southeast, and is slightly within the core 
radius ro of the ID surface brightness profile. 

The temperature of this feature is uncertain due to the 
low number of counts. The temperature from the adap- 
tive binning scheme is between kT ~ 3.5 — 5 keV, whereas 
the deprojected temperature must be kT > 5.5 keV, at 
least 2 — 3 times greater than the surrounding medium. 
We now consider whether this feature is consistent with 
a shock. 

The Rankine-Hugoniot shock-jump conditions are 



(7 + 



Pl__ 

Pi V2 (7 - l)Aff + 2 
Pi _ P2kT2/m _ 2-1 Ml - (7 - 1) 
P2 pikTi/m 7 + 1 

I2 ^ [(7 - l)Ml + 2] [27M2 - (7 - 1)] 



(7- 



(8) 
(9) 
(10) 
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where 7 = 5/3 is the adiabatic constant and Mi = v/cg is 
the Mach number. If we take T2/T1 w 2, then Mi w 1.9. 
Incorporating uncertainties from the thermal models and 
binning scheme, we find 

Afi-1.8-2.7. (11) 

Assuming Mi ^ 2, we would expect a downstream den- 
sity jump of p2 = 2.3pi, with a corresponding down- 
stream surface brightness ~5 times higher. However, 
farther downstream compressed material will relax, and 
the shock front is thin relative to the atmosphere. Thus, 
while we would expect to see a surface brightness discon- 
tinuity, it might be substantially smaller. 

In fact, we do not see such a discontinuity. The ID 
radial profile shows a hint of a surface brightness jump 
with S1/E2 ^ 1.4 in a region nearly cospatial with the 
temperature jump (Figure |4|), but there are insufficient 
counts to conclusively identify a discontinuity or detect 
the expected lag between the discontinuity and rise in 
temperature. Another hint is the elliptical surface bright- 
ness enhancement co-aligned with the jet axis (Figure [5]). 
This enhancement was earlier interpreted by us as the 
pre-outburst structure of the ICM, but it may result 
from the outburst. Unlike cold fronts in which a surface 
brightness discontinuity occurs in pressure equilibrium, 
the X-ray bright gas in the elliptical region has a some- 
what higher pressure than outside. A deeper observation 
is required to confirm the discontinuity. 

Another objection to the shock hypothesis is that the 
front is observed only on one side, but this can be ex- 
plained by the fact that the cluster atmosphere bright- 
ens to the northwest and tapers off to the southeast 
(Figure [2]). Shock-compressed gas will be most visible 
against a dim background, and the high-temperature 
knots near the northwestern lobe suggest that there is 
indeed shocked gas in the vicinity of the radio galaxy on 
the northern side. If the cluster atmosphere is denser in 
this direction, the shock wave will sweep up more ma- 
terial and advance more slowly per unit time than the 
southeastern front. The northwestern lobe is curtailed 
relative to the southeastern one, but ^^4 times more pho- 
tons are required to measure a density gradient on the 
relevant scales. 

If the high-temperature front is not a shock from the 
radio galaxy, its origin is not obvious. The radiative 
(free-free) cooling time of the hot gas in the region is 
more than 2 Gyr, but conduction due to hot electrons 
would also tend to spread out the thermal energy. It is 
possible this feature is a shock associated with a prior 
outburst that has since swept up material and become 
a large sound wave, but if it is traveling at the sound 
speed its distance from the active lobes would imply but 
a short delay between successive outbursts. It seems less 
likely that magnetic structures could produce a region 
(via anisotropic conduction) so much hotter than the sur- 
roundings. Without better data, we tentatively conclude 
that the feature is best explained as a shock associated 
with the current outburst. 

Active Lobe Age and Jet Power — Assuming that the ac- 
tive lobes represent a single outburst and inflate bubbles 
as large as the region of radio emission, we can make a 
rough estimate of the jet power from the size of the active 
lobes using standard arguments: the average jet power 



is the cavity energy divided by the age of the outburst. 

A note of caution is required because no cavity has 
actually been detected associated with the active lobes. 
Although it seems clear that radio jets must form bubbles 
in the ICM even when no cavity is detected, the filling 
factor of older cavities by ambient material is unknown. 
Assuming cylindrical cavities the volume depends only 
linearly on the projected length, so we sketch out order- 
of-magnitude estimates to place limits on the jet power. 

The energy required to blow a bubble in the ICM is 
Ecuv = lpV/{'~i — 1) = 4py, where p is the pressure, V 
the volume, and 7 = 4/3 the adiabatic constant for rela- 
tivistic electrons. The cavity enthalpy can be expressed 
as 

£^cav = 4py 47rrfob^/liobo"-gasfcrgas (12) 

where the gas subscript refers to the ICM. riobo ^ 20 kpc 
and /iiobo 75 kpc are estimates based on visual inspec- 
tion of the MOl maps and subject to projection and res- 
olution uncertainties; lower frequency maps suggest that 
a cylinder is a good representation. fcTgas 2.5 keV is 
obtained from 1-T fits to the region, and rigas is like- 
wise obtained from the spectral fit via the model emis- 
sion measu re. For assum ed geometry we follow the 
method in iMahdavi et al.l (2005), in which the emis- 
sion is assumed to come from a spherical shell with 
outer radius Ri and inner radius i?2- The line-of-sight 
depth is then 2{Ri — i?2)^/^. Near the active lobes, 
p ~ 10~^^ dyne cm~^. Hence, i?cav ~ 10^^ erg. Note 
that /iiobc is less certain than Hobe because of projection 
effects. Hence, if the wings were formed in the same 
outburst, the energy is only greater by a factor of ~2. 

A plausible range for the jet power can be obtained in 
view of the putative shock front. The presence of a shock 
suggests supersonic expansion, at least in the past — the 
radio jet is now more akin to that of an FR I source 
where subsonic expansion is expected. The Mach 
shock implies a source age of less than 75 Myr, but since 
shocks sweep up material and slow down, the initial ex- 
pansion was likely substantially faster. The well defined 
lobes that still give the radio galaxy an edge-brightened 
morphology suggest that the transition occurred recently, 
which is consistent with estimates of FR II and quasa r 
hfetimes of tens of Myr (jBird et al.ll2008t lMartinil[200l . 
The shock front is also still relatively close to the radio 
galaxy, so a scenario in which the jets are simply illumi- 
nating a very old source is disfavored. 

An average jet power of Pjot ^ 10^'* erg s^^, corre- 
sponding to a lobe age of 32 Myr, is consistent with 
the data. For reference, if NGC 326 had jet power 
commensurate with Cygnus A (Pjot ~ 10^^ erg s~-^; 
iStockton et al.|[l99^ . its lobes would be inflated in just 
0.3 Myr! There is no obvious cocoon and the source 
lacks hot spots, so this is clearly too high. On the 
lower end, transonic lobe expansion gives a jet power 
of Pjot ~ 5 X 10"*"^ erg s~^, which is inconsistent with 
the shock front. Pjot ^ 10'*'* erg s^^ is therefore prob- 
ably accurate to within a factor of 2 — 3, giving an age 
between 10-50 My r. This power is akin to Perseus A 
(jFabian et al.l[2006l ). 

4. DISCUSSION 
4.1. Reconstructing the History of NGC 326 
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In view of the high resolution X-ray data, we re- 
examine the formation history of NGC 326 and its long 
wings. The salient results of the X-ray analysis follow: 
(i) The very long wings in a relatively cool environment 
require a single outburst to have an unrealistic age if they 
expanded transonically or slower, (ii) The putative shock 
and high-temperature rims appear to be causally con- 
nected to the radio outburst, but are inconsistent with 
the present-day jet, which strongly resembles an FR I 
source: the jet decelerates and breaks up before reaching 
the ends of the lobes, and no terminal hot spots are ob- 
served. Thus, we do not expect that the strong backflows 
associated with the hot spots of FR II sources are present 
in NGC 326. Still, it is notable that the active lobes are 
well defined and bright, which led to the original classi- 
fication of the source as an FR II. (iii) The radio galaxy 
is situated at the region of brightest X-ray emission in 
the cluster, but the cluster emission shows structure on 
scales slightly smaller than the chip. The wings appear 
to follow an edge in this emission (Figure [2]). 

The data suggest that the outburst that produced 
the radio galaxy was much more powerful in the past. 
NGC 326 is very weak for an FR II source rel ative to 
its host galaxy luminosity (jCheung et al.l [2009| ) : if the 
present-day kinetic luminosity of the jet is in equiparti- 
tion with the radiative luminosity, the jet power is on 
the order of Pjet ~ 10''^ erg s~^. Although equipartition 
need not be the case, the jet would need to be extremely 
radiatively inefficient in order to drive the shock and in- 
flate the cavities if it is. The brightness of the jet near 
the core suggests this is not the case. If NGC 326 is the 
result of a powerful outburst that decayed, the shock and 
well defined lobes are natural consequences of the more 
powerful stage. 

In fact, the parameters of such an outburst fall within 
the accepted parameters for radio-galaxy power and life- 
time as well as the energy required in this system. Sup- 
pose that the jet power was initially Pjet ~ 10''^ erg s~^, 
more powerful than Perseus A but less powerful than 
Cygnus A, and decayed exponentially on an e-folding 
time of 10 Myr (thought to be a typical FR II lifetime; 
iBird et aLl i2008). Then the total energy input by the jet 
into the ICM is about 3 x 10^® erg, of which about two- 
thirds is injected during the FR II phase (0 — 10 Myr). 
After an e-folding time, the power of this hypothetical 
jet is only a few times higher than the observed radio 
power of the core, and the energy injected into the ICM 
is within a factor of 2 of our cavity enthalpy estimate, 
shock energy and cosmic ray heating notwithstanding. 
Hence, a fast-rise exponential-decay scenario seems to fit 
NGC 326. 

We can go a step farther since we know from the 
synchrotron cooling time that the plasma radiating 
at 5 GHz in the wings cannot be older than about 
40 Myr, and froni spec tral aging arguments (following 
IMvers fc Spangle^fTosI using the spectral data of MOl) 
that even the plasma radiating at 1.4 GHz is only about 
10 Myr old. If this plasma was deposited by flows 
from the primary lobes to the wings, the wing length 
of >100 kpc implies flow speeds of 3 — 10 x 10'^ km s~^, 
suggesting strong backflows (as in W95). Such back- 
flows are associated with the hot spots of FR II sources 
(e.g. iLeahv fc Williamslll984l : lAntonuccio-Delogu k SilS 



120101) . The radio galaxy must then have been strong 
enough to form hot spots no more than a few tens of Myr 
ago. Note that this argument does not rely on a specific 
mechanism of wing formation as long as the wings are il- 
luminated by active backfiows. Such an age is consistent 
with the location and tentative strength of the shock, 
which imply an outburst age of (significantly) less than 
75 Myr. If so, the length of the active lobes implies su- 
personic expansion, which is what we would expect from 
a jet driving a shock. 

But where did the wings come from? There are 
two general ideas in the literature. First, it is pos- 
sible that the wings are radio lobes which happen to 
grow in a different direction from the jet axis due to 
anisotropic pressure gradients in the hot atm ospheres 
(e.g. W95, Leahy & Williams igsl lCapettTet al. 2002). 
The other possibility is that the wings represent fossil 
lobes from a prior outburst in a different direction (e.g. 
iDennett-Thorpe et al.ll2Q02l: iMerritt fc Ekersl[200^ . As 
NGC 326 has been used as the poster-child of several 
ideas presented in the literature, we will evaluate both 
hydrodynamic and reorientation models. 

Powerful radio galaxies generate strong backfiows 
at the hot spot zones of the jet head. These 
flows are quite rapid (up to 20% of the jet speed; 
lAntonuccio-Delogu fc Silkri2010D despite the relatively 
slow advance of the jet head itself. W95 argue that these 
flows are buoyant in the ICM and, though rapid, may be 
deflected by a region of higher pressure. They find that 
flows of ^3000 km s~^ would be buoyant in the ICM of 
NGC 326 and propose that they formed the wings. How- 
ever, it is not likely that wings formed in this way ad- 
vance at a speed close to the backfiow speed. Simulations 
indicate that lobes are turbulent and that rapid back- 
flows experience shocks with in the lobe itself, dissipat- 
ing b ulk kinetic energy (e.g. iHodges-Kluck fc Reynoldd 
120111) . Although magnetic fields may collimate flows, 
mature magnetized lobes st ill become quite turbulent 
(|Huarte-EsDinosa et al.|[2(311l) . In our recent work specif- 
ically investigating the action of pressure gradients on 
radio lobes in pure hydrodynamic simulations, we found 
that wings expand at most transonically due to this dissi- 
pation of backfiow thrust. Hence, the presence of strong 
backfiows does not guarantee rapid wing expansion (in- 
deed, the internal sound speed of the lobes may be com- 
parable to the backfiow speed). In fact, it would be sur- 
prising if the simulations showed otherwise, since this 
would lead to wings significantly wider than the active 
lobes that also propagate faster than the jet head. 

This is a significant problem because the wings in 
NGC 326 are at least 100 kpc long; some XRGs have 
even longer wings and are in cooler atmospheres. Any ex- 
pansion mechanism that depends on the sound speed re- 
quires very old wings. Figure [8] shows the ages implied by 
(projected) wing lengths as a function of average expan- 
sion velocity for several XRGs with detected IGM/ICM 
emission. For reference, it would take ^200 Myr to in- 
flate the wings of NGC 326 transonically and 250 Myr to 
do so buoyantly. 

Another possibility is that the radio galaxy is born in 
a dense elliptical environment with the jet aligned with 
the major axis of the atmosphere. The cocoon formed by 
the spent jet material may then become overpressured 
if the jet advance is stalled by the dense atmosphere. 
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and at some point it will rupture along th e minor axis, 
thereby forming wings (jCapetti et al.ll2002l ). Before pres- 
sure equilibrium is established, the outflow will be super- 
sonic; afterwards, it will be akin to the buoyant backflow 
scenario. NGC 326 is not currently overpressurcd relative 
to its surroundings (MOl, .Worrall fc Birkinsh aw 20Q^, 
so despite the local IC M conforming r oughl y to the ge- 
ometry required in the iCapetti et alj ()2002D model, ex- 
pansion would currently be driven by buoyancy and bulk 
flows. As described above, this expansion would be sub- 
sonic. 

The circumstantial evidence therefore favors a reori- 
entation scheme. Detection of a bow shock ahead of the 
wings would be a strong indictment of the backflow mod- 
els, but there are insufhcient counts in the CXO image 
to perform an adequate search. We note that the tem- 
perature map in Figure [6] reveals bins to the northeast 
and northwest of the cavity that are 1.5 — 2 times hotter 
than the surrounding bins, and extraction of (low signal) 
spectra from bins half the size (~700 counts) suggests the 
ICM is indeed hotter closer to the cavity edges. About 
four times as many counts are required to determine if 
this is shock heating. 

On the other hand, the proble m with reorientation 
models, as noted by W95 and iGopal-Krishna et al.l 
(|2003f) . is that the wings have Z-shaped morphology and 
the eastern wing has an off-axis extension to the north- 
west (Figure [B- In respo nse to the Z-shaped wings, 
iGopal-Kr ishna et al.l (|2003| ) proposed a hybrid hydrody- 
namic and reorientation model in which a stream of ISM 
from a merging galaxy deflects the jet, thereby form- 
ing a ZRG. If a SMBH binary is formed by this merger, 
upon coalescence the jet axis will change and the ZRG 
becomes a fossil; the new lobes give the overall radio 
emission X-shaped mo rphol ogy. T his mod e l was pro- 
posed for NGC 326 by ICo j al-Krish na et al.l (|200l and 
updated for the same bvTZier (.2005^ with the sugges- 
tion that a stream of dense gas bent the jets at a radius 
between 25-50 kpc from the core. For reference, the 
projected separation between the nuclei of the dumbbell 
components in NGC 326 is ^6.4 kpc. 

A strength of the IGopal-Krishna et al.l ()2003f ) model 
is that it naturally explains the tendency for XRGs to 
fall near the FR I/II break (wit h the most recent con- 
firmation in iCheung et al.l [2OO90 . This is because de- 
flecting a jet requires (from purely hydrodynamic con- 
siderations) the pressure exerted by the jet on the ISM 
to be balanced by the pressure of the gas. For rea- 
sonable estimates of the density and bulk motion of 
streams of merging gas, highly relativistic jets cannot be 
stopped. Indeed, the non-winged FR II sources whose 
jets are deflected by interaction with dense gas also tend 
to have radio lu minosity close to th e FR I/II divide (e.g. 
PKS 2153-69; lYoung et all 120051) or hvbrid lobe mor- 
phology (e.g. 3C 321: Evans et al.ll2008l ). 

However, a jet weak enough to be deflected by the ram 
pressure of the ISM would likely inflate lobes at some 
fraction or low multiple of the sound speed (as is thought 
to be the case with FR I sources). Tying expansion to 
the sound speed suggests very long radio galaxy lifetimes 
(Figure [S]) . Apart from fueling the AGN steadily, this 
timescale would require the bending agent to be in con- 
tact with the jet at favorable orientations on both sides of 
the galaxy for up to hundreds of Myr. A gas stream 25- 
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Fig. 8. — Wing expansion times as a function of the Mach number 
for several XRGs with X-ray data. Note that because galaxy group 
atmospheres are relatively cool, modest expansion rates may still 
be very supersonic. Wing formation models that grow wings at a 
fraction or multiple of the sound speed require old wings. 



50 kpc from the AGN with bulk motion of ~ 2G0 km s~^ 
(comparable to streams in similar systems; iZieii |2005[) 
would travel perhaps tens of kpc in this time, meaning 
the stream must be very long and fortunately oriented 
in order to keep the jet bent at roughly the same angle 
at the same radius from the core. Although conceivable, 
this is a very special circumstance. Another objection in 
the same vein is that the time between the production 
of these streams of rotating gas and the coalescence of a 
SMBH binary that engenders the present X-shaped mor- 
phology must also be up to hundreds of Myr or the wings 
would not have gotten so long, regardless of the stream 
geometry. The time required for a SMBH binary to coa- 
lesce is unknown, but one might imagine that if the radio 
galaxy were driven by a component of the SMBH binary 
over 100 Myr, accretion from the same disk of gas would 
cause the spins to become aligned prior to the merger and 
forest all a spin-flip upon co alescence (iBogdanovic et al.l 
I2007f ). It is possible that the IGopal-Krishna et al.l (|2003D 
model is essentially correct, but it faces two hurdles: 
there is presently no evidence for the required streams 
in NGC 326 and the streams, if found, must be able to 
bend a powerful jet. 

It is possible that the Z-shaped morphology formed 
after reorientation. The buoyant backflow and over- 
pressured cocoon models both take advantage of the 
anisotropy of the atmosphere to push the radio lobes in 
the right direction. This motion would still occur with 
fossil relics. Suppose that, while active, the lobes of the 
wings formed an axis across the nucleus, but since the ex- 
tinction of the jet in that direction evolve buoyantly. If 
the major axis of the pocket of ICM they inhabit is some- 
what off-axis from the major axis of the radio galaxy (but 
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still close enough to satisfy the lCapetti et al.l 120021 cor- 
relation) , buoyancy will tend to lift the wings away from 
the higher pressure core in opposite directions. The cen- 
ters of the wings are about 25 kpc away from a line that 
would connect them to the nucleus, implying a rise time 
of less than 50 Myr. This is consistent with the age of the 
active lobes. There is evidence that the wings do respond 
to the pressures in the ICM: both appear to be follow- 
ing a ridge of brighter ICM emission (Figure [5]), forming 
a large but shallow arc. Assuming the wings were bent 
into this shape after the jet ceased, the radius of curva- 
ture compared to a straight line through the wings sug- 
gests buoyancy would produce this morphology within 
40-70 Myr depending on the distance from the core (this 
estimate does not take into account spatial variations 
in the buoyant terminal velocity due to an anisotropic 
medium). However, the plume to the northwest of the 
eastern wing (Figure [IJ would take longer to rise and is 
not well explained by this scenario. 

4.2. Wings and Radio-Mode Feedback 

In any plausible wing formation scenario, the very long 
wings of NGC 326 are presently evolving passively in the 
ICM, i.e., as buoyant structures subject to hydrodynamic 
forces. They are there fore akin to "dead" radio galaxies 
(|Revnolds et al.1 12002D in which t he jet has ceased (in 
the model o f lMerritt fc E kCTi'200^. this is precisely what 
they are). Dead radio ga laxies arc notoriously rare (for 
one example, see lWorrall et al.l l2007). but crucial to the 
paradigm of "radio-mode" AGN feedback. 

In radio-mode feedback, radiative cooling in the ICM 
is regulated by the condensation of gas onto the central 
galaxy, fueling an AGN that drives energetic outflows 
into the ICM, heating it. X-ray cavities associated with 
radio lobes provide incontrovertible evidence that these 
outflows exist, interact with the ICM by inflating bubbles 
of relativistic plasma in it, and deposit suffic ient energy 
to re gulate cooling (for a recent review, see iGitti et al.l 
l2011f) . However, it is not clear that these outflows do 
regulate cooling, since we do not know how this energy 
is thermalized and isotropized or on what timescales this 
occurs. 

Several heating mechanisms have been identified. 
These include active heating through strong or weak 
shocks; adiabatic expansion of a jet-blown cavity; and 
damped sound waves associated with the cavity infla- 
tion, and passive heating such as entrainment of mate- 
rial in a buoyant pl ume or infall of displaced g as behind 
a buoyant bubble ijMcNamara fc NulsenI 120071 and ref- 
erences therein). Although X-ray observations demon- 
strate that all of these occur, their relative importance 
(which may be a function of time, environment, and jet 
power) is unknown. A key problem is that the life cy- 
cle of AGNs (which itself may be a function of redshift 
and whether the host galaxy is in a group or a cluster) 
is not understood. For example, an AGN with a short 
duty cycle punctuated by brief, powerful outbursts may 
be less effective at heating than an AGN that more or 
less continually blow s small bubbles ( "effervescent" heat- 
ing; iBegelmanl 120011 ) . Dead radio galaxies are therefore 
of interest, since their frequency is a measure of the AGN 
duty cycle and the rate at which they heat the ICM is a 
measure of the importance of passive heating. 

Thus, wings may be of interest as probes of feedback. If 



XRG wings are dynamically similar to dead radio galax- 
ies, they offer a few additional advantages as laboratories 
of dead sources. First, dead radio galaxies are most fre- 
quently detected as "ghost" cavit i es un a ssociated with 
GHz emission (e.g. F abian et all 120001 : iRaffertv et all 
[200 6; Dong et al. 2010), but ghost cavities are typically 
detected near the group or cluster core where the surface 
brightness is high. It can therefore be difficult to dis- 
tinguish between features associated with old outbursts 
and the active ones. Wings, on the other hand, can ex- 
tend more than 100 kpc from the AGN and in a very 
different direction. Second, it is difficult to age ghost 
cavities, whereas a better understanding of wing forma- 
tion would yield the ages of both the active lobes and 
wings in XRGs. Third, winged sources are typically lo- 
cated in group rather than cluster environments. Until 
recently, feedback in groups has been largely neglected, 
but non-gravitational processes are even more important 
in these systems than in clusters, and the total number of 
galaxies in groups far outweighs the cluster population. 

Our examination of the formation history of NGC 326 
suggests another advantage to XRG wings: if the rea- 
son XRGs are typically found near the FR I/II break 
because they represent decaying outbursts, the total ac- 
tive and passive energy injection rates over the history 
of the source can be measured in the same system. The 
deep X-ray observations required for this kind of mea- 
surement only exist in a handful of clusters with both 
active lobes and ghost cavities (the most prominent, of 
course, is Perseus: , Fabian et al.112006, '). Since groups are 
much fainter, it is essential to have good targets. 

The CXO observation of NGC 326 is not deep enough 
to carry out heating measurements, but the surface 
brightness of features detected suggests that an exposure 
with four times the counts (twice the signal) would suf- 
fice. For example, the high-temperature front has about 
4,000 counts in its vicinity. To determine if it is a shock 
and measure its strength to 10% accuracy requires sur- 
face brightness values and deprojected temperatures to 
at least this accuracy in semi-annular bins about 15 arc- 
sec wide at a radius of 50 kpc from the nucleus. This 
would require about 15,000 counts per bin. Since strong 
shocks are much more efficient at irreversibly heating the 
ICM than weak ones, the shock strength is of interest. 

Likewise, the cavity in the wing would be detected at 
the 7(7 level in an exposure of this duration, and cav- 
ities associated with the other lobes (which we assume 
to be there) would be easily detected, barring patholog- 
ical projection effects. In the wing, this quality of de- 
tection would place tight constraints on the cavity size, 
important for determining the cavity stability and plac- 
ing boundaries from which to measure the amount of 
entrained material, uplifted material, and temperature 
or metallicity gradients. An image with twice the signal 
would allow for adapative binning at scales significantly 
smaller than the lobe size with temperature measure- 
ments accurate to 20%. Specific entropy (a surer mea- 
sure of irreversible heating than temperature) could also 
be measured on these scales. Metallicity could be deter- 
mined in slightly larger bins that would still be smaller 
than the lobe size, enabling a measurement of metallicity 
gradients around and behind the cavity. Given an ac- 
curate measurement of the terminal buoyant vel ocity of 
the bubble, the displacement heating rate (|Birzan et al.l 
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I2004D could then be measured. 

Thus, while the extant CXO image lacks the depth 
to take advantage of the measurements sketched here, 
it demonstrates that such measurements are well within 
reach. The utility of XRG wings depends on high-quality 
X-ray data and better understanding wing origin, so flux 
estimates in other systems are important. Only 13 XRGs 
have CXO or XMM-Newton data, and only ^8 of these 
have significant I GM/ICM emission (mode l flux esti- 
mates provided in IHodges-Kluck et al.|[20 Tot). Only two 
of these have moderately deep exposures (NGC 326 and 
4C -1-00.58). The CXO image presented here suggests 
that for larger XRGs, the higher effective area of XMM- 
Newton makes it a better instrument for these types of 
measurements . 



5. SUMMARY AND CONCLUSIONS 

We report on a 100 ks CXO image of NGC 326, an 
XRG with wings 100 kpc long that resides in a poor 
cluster. Owing to its dramatic wings, NGC 326 has fre- 
quently been the su bject of studies atte r apting to explain 



wing origin fW95: iWirth et al 



2002; Gopal-Krish na et all 120031: iDennett-Thorpe et al 



[1981 iMerritt fc Ekers 



2002t iZien i2005c iCaoetti et all I2002t ILal fc Raol 120051) 



The CXO image provides an important step forward in 
evaluating these hypotheses due to the identification of 
several X-ray features associated with the radio galaxy. 
These include a high-temperature front to the southeast 
of the radio galaxy (Figure [6]), which we assume to be a 
shock front, high-temperature knots of emission around 
the radio galaxy akin to cavity rims (Figure [5]) , and a 
cavity associated with the eastern wing (Figure [T]) • The 
high angular resolution of CXO compared to ROSAT 
also allows us to remove point sources and character- 
ize the broader ICM, revealing an X-ray peak associated 
with the northern lobe of the radio galaxy rather than 
the ISM and an elliptical enhancement around the radio 
galaxy (first identified in iHodges-Kluck et al.l[20Tol ). 

By characterizing the immediate environment of the 
radio galaxy and taking into account spectral aging of the 
radio emission in the wings, we can construct a reason- 
able history for the AGN outbursts. The jets, lack of hot 
spots, and low radio power (lO''^ erg s"^) indicate a weak 
FR I- like source today, but the putative shock front, lobe 
size and morphology, and GHz emission at the wing tips 
indicate a powerful FR II source no more than ^30 Myr 
ago. We therefore suggest the kinetic jet power could be 
described as Pjct — W^^e 



-tjta 



erg s~^, where t^ — Wi Myr 
for the current outburst. This function would account for 
the observed features as well as produce about the right 
amount of cavity enthalpy. The energy required to infiate 
cavities the size of the lobes places the current kinetic jet 
power between 5 x 10**^ — 5 x 10"*^ erg s~^. The radio 
power suggests that the flow is therefore still radiatively 
inefficient. 



NGC 326 is the best studied XRG, and the variety of 
wing formation models associated with it is a testament 
to the difficulty of proving any one scenario. The sound 
speed in the ICM is only 500 km s"^ (0.5 kpc Myr'^), 
implying that the long (100 kpc) wings must have ex- 
panded supersonically. This disfavors the backflow mod- 
els, which have yet to convincingly demonstrate that in- 
ternal lobe flows can generate sustained supersonic ex- 
pansion; simulations indicate that they become turbu- 
lent and dissipate their thrust before reaching the wing 
edges. In fact, any model in which the wing growth is 
related to the sound speed is disfavored in systems with 
long wings (Figure [S]). The most likely alternative is that 
the wings are fossil relics illuminated by the present out- 
burst, but Gopal-Krishna et al.l ()2003D note that a simple 
spin-flip would not explain the Z-shaped morphology of 
the wings. However, their model relies on bending a mod- 
erately weak jet, which would not be expected to propel 
hypersonic lobes. As fossil radio lobes would be subject 
to hydrodynamic forces, it is possible that the Z-shaped 
morphology of the wings is produced by buoyant evo- 
lution in an anisotropic atmosphere. It would take less 
than 50 Myr for the wings to rise buoyantly away from 
an axis connecting them through the nucleus to their 
present-day location, which is consistent with the age of 
the active lobes and the age required for the wings to re- 
spond to the anisotropy of the large-scale cluster emission 
to form the understated C-shaped arc (Figure [5|). How- 
ever, this scenario fails to explain the plume on the east- 
ern wing. It is also notable that the tendency for XRGs 
to have j ets co-aligned with the major axes of their host 
galax ies (jCapetti et al.|[20Q2t iSaripalli fc SubrahmanvanI 
|2009() and hot atmospheres (Hodgcs-Kluck et al. 2010) 
extends to XRGs with very long wings. This is unex- 
plained by any reorientation model, and remains com- 
pelling evidence that the environments of most XRGs 
help determine their morphology. 

XRGs provide an important target for characterizing 
AGN feedback through the study of "ghost" cavities, 
both as pathfinders (they illuminate the cavities in the 
GHz bands) as well as targets which are readily distin- 
guished from the active outburst. The heating rate as- 
sociated with dead outbursts can then be measured in 
isolation, and compared to the timescales required to 
thermalize and isotropize the energy in the ICM. In cases 
such as NGC 326, where we believe most of the energy 
associated with the active outburst has already been de- 
posited, we can use the AGN lifetime to estimate an av- 
erage heating rate, assuming some duty cycle. XRGs are 
therefore interesting beyond the origin of their bizarre 
morphology. 

The authors thank the anonymous referee for helpful 
comments. GSR thanks support from the National Sci- 
ence Foundation under grant AST0908212. 
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